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FOREWORD 

This  in-house  effort  was  accompl i shed  under  Work  Unit  24041QA1,  entitled 
“Aerodynamic  Design  and  Analysis  Methods."  The  work  reported  herein  was 
performed  during  t.hp  period  March  1906  to  November  1986  under  the  direction  of 
the  author,  James  R.  Sirbaugh  (AFWAL/FIMM) . 
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SECTION  I 
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INTRODUCTION 

The  solution  of  the  F.uler  equations  for  aircraft  flow  fields  involves  two 
major  problems:  grid  generation  and  flow  equation  solution.  A  grid  must  be 
generated  for  each  new  configuration  t.o  be  studied.  The  grid  must  accurately 
model  the  configuration  surface  geometry  and  provide  sufficient  grid 
resolution  in  t-h«  reg'i on/around  the  conf iuuration  to  capture  the  flow  details. 
Due  to  current  computer  speed  and  memory  limits,  a  computational  aerodynamics 
cnrle  must  be  limited  in  the  number  of  grid  points  it  uses,  and  therefore  the 
qric  point,  locations  must  be  chosen  carefully.  The  following  pages  describe 
the  development  and  application  of  a  delta  wing  grid  topology. 

Grid  generation  methods  should  be  tailored  according  to  both  the  physics 
of  the  flow  and  the  flow  equation  solution  method.  The  primary  flow  charac¬ 
teristic  of  the  delta  wing  is  the  leading-edge  separation  that  rolls  up  into  a 
vortex.  The  vortex  position,  size  and  strength  are  dependent  on  many  factors, 
one  of  which  is  the  leading-edge  shape.  A  sharp  leading  edge  will  generally 
produce  a  vortex  originating  at  the  edge,  while  a  round  leading  edge  may 
produce  a  separation  point  more  inboard  on  the  upper  surface  thus  increasing 
the  strength  of  the  vortex.  The  leading-edge  grid  must  be  fine  enough  to 
permit  the  flow  solver  to  capture  the  flow  gradients  contributing  to  leading- 
rdqe  separation  in  order  to  accurately  predict  vortex  core  position  and 
strength.  Downstream  of  the  wing,  the  vortex  core  will  diffuse  to  a  degree, 
and  may  even  burst  or  breakdown.  If  a  computational  method,  such  as  an  Euler 
solver,  is  going  to  predict  the  occurrence  of  vortex  burst,  the  downstean  grid 
must  adequately  resolve  the  vortex  flow.  Fvpn  for  cases  where  burst  is  not  a 
consideration ,  the  grid  immediately  downstream  of  the  wing  must  be  dense 
enough  to  accurately  predict  the  trai ling-edge  flow  field. 


The  Euler  equation  solver  used  in  this  work  is  FIP57  (Ref.  1).  While 
fl  057  has  previously  been  applied  to  delta  wing  configuration,  t.he  built-in 
grid  generation  in  FL057  is  inadequate  for  delta  wings.  The  grid  does  riot 
mode  1  the  wingtip  and  produces  a  downstream  grid  which  is  highly  skewed  lor 
configurations  with  large  leading-edge  sweeps.  A  few  aspects  of  FL057  should 
be  kept  in  mind  when  considering  a  grid  topology  to  be  coupled  to  the  Euler 
equation  solver.  A  finite-volume  and  central -difference  scheme  is  used  by  the 
method.  The  finite-volume  method  does  not  require  grid  lines  to  be  continuous 
in  slope,  or  orthogonal;  thus,  the  grid  can  be  generated  by  a  variety  of  ways, 
including  those  which  do  not  directly  enforce  orthogonality.  The  total  grid 
may  also  be  broken  up  into  several  regions  and  each  region  generated 
separately,  as  long  as  grid  lines  are  continuous  across  region  boundaries.  A 
highly  skewed  grid,  however,  will  generate  inaccurate  flow  field  solutions; 


i.d.,  a  fact  that  must  be  kept  in  mind  when  considering  nonorthogonal  grids. 

The  central -di fterence  scheme  is  used  to  calculate  mass,  momentum,  and 
energy  fluxes  across  each  grid  cell  face  by  averaging  the  flow  variable  values 
adjacent  to  the  cell  face.  This  calculation  procedure  assumes  that  the 
distance  between  the  cell  face,  where  the  flux  is  calculated,  and  the  cell 
center,  where  the  flow  variables  are  calculated,  is  approximately  the  same  for 


neighbor i ng  cells.  Should  the  grid  point  spacing  vary  extremely  rapidly  or  in 
a  nonsmooth  fashion,  the  simple  averaging  ot  variables  will  produce  an  inaccu¬ 
rate  approximat ion  to  the  cell  face  flux.  The  far-field  grid  region  typically 
has  rapid  grid  stretching;  however,  the  flow  solution  accuracy  is  not  signifi¬ 
cantly  degraded  in  this  region  due  to  the  nearly  uniform  flow  at  the  far 
field. 

Boundary  conditions  are  applied  in  the  FL057  Euler  method  in  a  variety  of 
ways ,  some  of  which  require  an  extrapolation  of  interior  flow  variables  t.o  the 
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grid  boundaries.  Although  the  boundary  conditions  do  not  explicitly  require* 
orthogonal  grids,  it  is  reasonable,  whpn  considering  the  boundary  condition 
extrapolation,  that  orthogonality  at  the  grid  boundaries  would  improve  the 
code's  accuracy. 
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SECTION  II 


GRID  TOPOLOGY 

The  new  delta  wing  grid  developed  during  this  effort  is  based  on  a  previ¬ 
ously  developed  two-dimensional  ( 2-D )  grid  topology  for  equilateral  trianyles 
(Ref.  ?) ,  Fig.  I.  The  leading  edge  and  wingtip  are  considered  to  be  a  single 
side  of  the  equi 1 ateral -triangle  grid.  The  singularity  point  is  positioned  on 
the  root  chord  and  therefore  should  not  be  in  any  strong  flow  gradients  region 
for  most  problems  of  interest.  Figures  ,°a-c  show  the  wing  surface  topology 
tor  a  65°  leading-edge  sweep,  cropped  delta  wing.  This  topology  permits  a 
uniform  grid  point  spacing  and  cell  size  along  the  entire  length  of  the 
leading  edge.  Grid  cells  can  be  clustered  forward  on  +he  wing  by  moving  the 
singularity  point  forward,  thus  permitting  adequate  grid  resolution  to  capture 
the  strong  flow  gradients  in  the  apex  region. 

The  new  delta  wing  grid  is  intended  primarily  for  use  with  the  FL0b7  Euler 
analysis  code.  In  order  to  couple  the  new  grid  topology  to  the  existing  FL057 
Euler  code  with  as  few  changes  to  FL05/'  as  possible,  the  following  grid  index 
system  was  adopted.  The  grid  lines  that  wrap  across  the  wing  surface  and 
around  the  leading  edge  of  the  wing  will  be  referred  to  as  the  I-varying  (K 
constant)  grid  lines.  (See  3a.)  The  1=1  grid  points  are  located  on  the 
wing's  lower  surface  root  chord.  The  1  -  I  grid  points  are  located  on  the 
wmg's  upper  surface  root  chord.  The  line  forming  the  forward  part  of  the 
root  chord  starts  at  the  lower  surface  singularity,  varies  in  I  and  moves 
forward  around  the  leading  edge  and  proceeds  aft  to  the  upper  surface  singular¬ 
ity.  This  grid  line  is  varying  in  I  and  has  a  constant  Y  value  of  one.  The 
t -varying  family  of  curves  start  at  the  forward  plane  of  symmetry  with  a  value 
nf  one  and  follows  paths  aft  to  the  downstream  exit  plane.  Thus,  the  winq's 
trailing  edge  is  formed  by  a  grid  line  varying  in  I  and  a  constant  K  value. 
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To  form  a  three-dimensional  (3-D)  grid,  the  C-H  grid  topology  in  FL057 
was  partially  merged  to  the  delta  wing  grid  topology  to  form  a  wrapped  0-H 
grid.  Fach  0  wrap  followed  an  I-vurying  g»~i d  line,  with  the  0  wraps  being 
almost  spanwise  cuts  near  the  wing's  trail ing  edge.  A  typical  0  wrap  from  the 
?-[  grid  is  shown  in  Figs.  4a  and  b.  The  wing  surface  was  formed  by  the  J  -  1 
grid  points  and  the  J  grid  points  form  the  outer  shell  of  the  3-D  grid. 

The  delta  wing  grid  topology  is  not  totally  new  and  is  easy  to  follow  if  the 
viewer  considers  the  grid  to  be  a  series  of  2-D  0  grids  stacked  together  to 
form  the  3-D  grid.  Figures  3a  and  b  show  the  pi ane-ol -symmetry  grid  with 
met/  directions  indicated. 

To  couple  the  3-D  delta  wing  grid  to  FL057,  the  boundary  conditions  to  be 

applied  on  each  face  of  the  3-D  computational  cube  must  be  determined.  The  J 

=  1  side  is  the  wing  surface  for  all  K  values  less  than  or  equal  to  the  K 

value  of  the  wing's  trailing  edge.  For  J  =  1  and  K  values  greater  than  K 

trailing  edge  value,  the  grid  collapses  to  a  slit  line  and  a  flow  through 

boundary  condition  is  used.  The  J  =  Jmax  side  is  the  far-field  grid  and 

appropriate  in-flow/out-flow  boundary  conditions  apply.  The  plane-of-symmetry 

grid  is  formed  by  the  I  =  1,  I  =  I  and  thp  K  -  1  sides  of  the  computational 

grid.  The  only  remaining  side  of  the  computational  spacp  is  the  K  =  K  side 

max 

and  this  is  the  downstream  boundary  grid  where  appropriate  out-flow  boundary 
conditions  are  applied. 
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GRID  GENERATION 


The  delta  wing  grid  used  in  the  Euler  computation  was  developed  in  three 
phases:  2-D  boundary  grids,  3-D  interior  grid,  and  algebraic  downstream  grid. 
The  2-D  grids  are  the  six  sides  of  the  3-D  grid,  consisting  of  the  winy  upper 
and  lower  surfaces,  symmetry  plane  (3  sides),  far  field  outer  shell,  and  the 
trai  1  ing-edge  plane.  The  symmetry  plane  is  generated  as  a  single  surface,  hut. 
is  actually  three  sides  of  the  computational  space  as  described  in  the  grid 
topology  section.  It  is  desirable  to  directly  control  grid  spacing  and  shape 
upstream  of  the  root  chord,  and  spanwise  beyond  the  wingtip.  To  control  the 
wingtip  region,  the  delta  wing  grid  is  divided  by  a  2-D  surface  grid  along  the 
wing  trailing  edge  extending  outwards  to  the  far-field  grid.  Grid  spacing  on 
the  2-D  divider  grid  is  specified  in  both  spanwise  and  vertical  directions, 
thus,  providing  the  desired  grid  control  in  the  wingtip  region.  The  3-D  grid 
upstream  of  the  wing  trailing  edge  is  complex  in  topology  and  is  created  by 
the  3-D  elliptic  grid  generation  program,  which,  when  properly  applied, 
produces  a  smoothly  varying  grid.  The  downstream  grid,  on  the  other  hand,  is 
essentially  2-D  in  topology  and  is  created  by  an  algebraic  method.  The 
junction  between  the  two  subgrids  is  not  continuous  in  regards  to  grid  line 
slope,  but  a  mild  variation  in  slope  continuity  will  not  prooi.ce  significant 
errors  in  a  finite  volume  flow  solver. 

TWO-DIMENSIONAL  BOUNDARY  GRIDS  (2-D 


T he  wing  surface  topology  as  previously  described,  is  not  unique,  but  is 
very  demanding  of  the  grid  generation  method.  A  desirable  feature  of  such  a 
method  is  the  ability  to  specify  the  grid  spacing  along  the  root  chord, 
trailing  and  leading  edges,  then  assimilate  this  spacing  smoothly  into  the 
interior  of  the  grin.  A  method  for  generating  quasi  2-D  grids  on  curved 


surfaces  by  Thomas  (kef.  3)  exhibits  the  desired  feature.  The  quasi  2-D 
method  was  derived  from  the  3-D  elliptic  grid  generation  method  by  Thompson 
(Pet.  2),  by  rewriting  the  3-D  equation  for  grid  points  being  constrained  to  a 
3-D  surface.  Thus,  the  2-D  method  accounts  for  grid  point  spacing  and 
position  based  on  a  3-D  position  vector,  not  a  2-D  position  vector.  The  wing 
leading-edge  region  is  where  the  added  dimension  makes  the  biggest  difference. 
An  additional  tern:  in  the  quasi  2-D  method,  but  not  found  in  the  3-D 
equations,  draws  grid  points  toward  regions  of  high  curvature,  thus,  producing 
a  more  accurate  modeling  of  the  curvature.  The  grid  spacing  in  the  interior 
of  the  grid  is  controlled  by  the  grid  point  spacing  on  the  boundaries  by 
adding  source  terms  to  the  elliptic  equation.  The  source  term  values  are 
computed  solely  from  the  boundary  grid  point,  spacing  and  then  interpolated 
inlo  the  interior  grid.  The  final  product  of  this  method  is  a  grid  that  has 
smoothly  varying  grid  point  positions  and  reflects  the  grid  point  spacing  of 
the  boundaries.  Appendix  A  lists  the  quasi  2-D  elliptic  equations. 

MING  SURFACE 

The  wing  upper  and  lower  surfaces  are  generated  separately  and  then 
joined  together.  Grid  point  spacing  is  specified  at  four  points  on  the  wing 
surfaces,  as  shown  in  Fig.  2c  with  TI,  T2,  T3,  and  T4  as  input  parameters. 

The  root  chord  is  first  spline  fitted  from  fhe  input  airfoil  shape.  Grid 
points  are  then  distributed  along  the  root  chord  based  on  arc  length,  with  the 
leading-edge  spacing  being  Tl  percent  of  thp  total  root  chord  arc  length  and 
the  tro i 1 i ng-edge  spacing  being  T2  percent  of  the  total  arc  length.  A  hyper¬ 
bolic  tangent  stretching  routine,  developed  by  Vinokur  (Ref.  4),  is  used  to 
determine  the  remaining  grid  point  positions.  The  grid  point  spacing  along 
the  trailing  edge  is  controlled  by  the  T3  and  T4  parameters.  The  spacing  at 
the  wingtip  is  T3  times  Tl;  therefore,  a  value  of  1.0  for  T3  will  create  an 


constant,  grid  spacing  along  the  leading  edge  over  the  entire  length  of  the 
leading  edge  and  wingtip.  Even  spacing  is  used  along  the  combined  length  of 
the  leading  edge  plus  wingtip.  The  grid  spacing  along  the  trailing  edge  at 
the  root  chord  is  T4  percent  of  trai 1 ina-edge  arc  length.  The  remaining 
i rai 1 ing-edge  points  (between  T3  and  T4)  are  distributed  based  on  the  hyper¬ 
bolic  tangent  stretching  function. 

The  wing  surface  is  generated  using  the  quasi  2-D  method  and  the  boundary 
spacing  as  just  described.  The  thickness  of  the  wing  at.  each  grid  point  is 
linearly  interpolated  from  airfoil  cross  sections  provided  as  user  input.  The 
wingtip  is  rounded  over  the  last  2  percent  of  semispan  as  seen  in  Fig.  5. 

PLANE  OF  SYMMETRY 

The  plane-of-symmetry  grid  (or  wall  in  a  wing-on-a-wal 1  calculation)  is 
created  by  the  same  2-D  elliptic  method  as  the  wing  surface,  simply  by 
rotating  the  symmetry  plane  into  the  wing  plane  and  setting  all  thickness 
values  to  zero.  The  symmetry  plane  is  divided  into  upper  and  lower  halves, 
and  a  separate  grid  is  generated  for  each  half.  The  two  halves  are  rejoined 
and  the  joining  line  is  adjusted  to  give  smooth  grid  lines  across  the  joint. 
The  upper  and  lower  halves  of  the  symmetry  plane  with  the  sides  numbered  for 
reference  are  shown  in  Fig.  6a-d.  Side  two  is  the  wing  root  chord  with  grid 
point  spacing  defined  by  the  wing  surface  grids.  Sides  one  and  four  were 
straight  lines  with  grid  spacing  at  the  root  chord,  far  field  and  overall 
lengths  of  these  lines  specified  by  the  user  as  input  parameters  T5,  T6,  and 
T7  on  side  ore  and  T8,  T9,  and  T10  on  side  four.  The  remaining  grid  points 
along  sides  one  and  four  are  distributed  using  the  hyperbolic  tangent 
stretching  routine. 

Once  the  upper  and  lower  halves  of  the  symmetry  planes  are  created,  the 
grid  points  on  side  one  are  adjusted  to  obtain  continuous  second  derivatives 


B 
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in  grid  spacing  lor  the  family  of  lines  crossing  side  one.  By  using  two  grid 
points  on  one  side  of  side  one  and  one  grid  point  of  the  same  curve  on  the 
other  side,  second  derivatives  in  grid  spacing  can  be  computed  on  both  sides 
of  side  one.  By  averaging  the  two  values  obtained,  a  first-order  accurate 
approximation  of  the  second  derivative  in  grid  spacing  is  obtained  for  points 
on  side  one.  Using  this  approximation  and  the  position  of  the  grid  points 
immediately  on  either  side  of  side  one,  the  grid  points  on  side  one  are 
repositioned  to  satisly  the  second-order  derivative  equation  written  for  side 
one  points.  The  result  of  this  repositioning  is  a  continuous  second 
derivative  in  grid  spacing  across  the  interface  of  side  one. 

The  2-D  elliptic  method  used  does  not  directly  enforce  orthogonality  at 
the  boundaries.  The  Si  term  in  the  source  terms  influences  the  curvature  of  a 
grid  line  (Ref.  3  and  Appendix  A).  The  source  term  is  normally  computed 
automatically  from  the  boundary  grid  point  position  and  spacing.  By  locally 
adjusting  the  value  of  SI,  the  boundary  grids  are  made  more  nearly  orthogonal. 
The  SI  parameter  was  set  by  trial  and  error,  for  the  grids  included  in  this 
report  at  the  wing  leading  edge  in  the  plane-of-symmetry  grid,  and  at  the 
wingtip  in  the  trai 1 i ng-edge  grid  to  improve  orthogonal ity  at  thp  wing  surface 
beyond  that  which  is  obtained  from  using  the  automatic  controls. 

TRAIL ING-EDGE  GRID 

The  ?-D  trai 1 ing-edge  grid  divides  the  total  3-D  physical  grid  into  two 
parts:  the  portion  upstream  of  the  wing  trailing  edge  and  the  portion  down¬ 
stream  of  the  trailing  edge.  This  grid  is  created  in  essentially  the  same  way 
as  the  symmetry  plane  grid,  but  the  trail ing-edge  shape  and  thickness  is 
modeled  instead  of  the  root  chord.  The  upper  and  lower  halves  of  the  grid  are 
created  separately.  The  halves  are  joined  together  and  the  joint  line  is  then 
smoothed.  Indicated  in  Fig.  7  is  the  trai 1 ing-edge  grid  side  numbering  and 
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spacing  parameters  used  in  the  computer  code.  Note  the  use  of  T8  and  T9  to 
control  spanwise  grid  spacing. 


Side  two  is  the  wing's  trailing  edge  and  is  generated  during  the  wing 
surface  grid  calculations.  Sides  one  and  four  are  straight  lines  with  end 
point  grid  spacing  and  overall  lengths  specified  as  input  parameters.  The 
remaining  side  one  and  side  four  grid  points  are  distributed  by  the  1-D  hyper¬ 
bolic  tangent  stretching  routine.  Side  three  is  elliptic  in  shape  and  has 
even  grid  point  spacing.  Figures  8a  and  b  show  a  typical  trai 1 ing-edge  grid. 
Even  though  this  grid  is  generated  in  2-D,  the  wing  trailing  edge  need  not  be 
a  straight  line. 

FAR-FIELD  CAP  GRID 

The  last  grid  required  to  complete  the  boundaries  of  the  3-D  physical 
space  is  the  far-field  cap  grid.  This  grid  is  an  ellipsoid  that  is  created 
using  the  quasi  2-D  elliptic  method.  The  far-field  cap  is  identical  to  the 
wing  surface  in  topology  and  side  numbering  (Figs.  9a  and  b).  All  sides 
except  the  elliptical  shaped  side  three  are  defined  by  sides  of  the  other 
quasi  2-D  surfaces.  An  even  grid  spacing  for  the  entire  farfield  cap  is  used 
on  side  three. 

THREE-DIMENSIONAL  GRID  (3-D) 

The  3-D  grid  method  used  to  generate  the  delta  wing  grid  was  developed  by 
Thomas  (Ref.  3).  This  method  is  also  the  basis  for  the  2-D  method  previously 
described.  The  standard  3-D  elliptic  equation  was  extended  to  Poisson's 
equation  by  inclusion  of  source  terms  that  caused  a  mimicking  of  the  2-D 
boundary  grid  point  spacing  in  the  interior  of  the  3-D  grid.  Appendix  B 
outlines  the  3-D  elliptic  equations  used  in  this  method.  Since  all  grid 
spacing  control  is  evaluated  from  the  2-D  boundary  grids  in  this  method,  the 
3-D  grid  generator  can  be  applied  to  a  variety  of  topologies  without  any 


special  considerations.  The  grid  singularity  found  in  this  delta  wing  grid 
topology  did  prevent  standard  calculation  of  the  SI  curvature  term  at  the 
singularity.  The  SI  term  is  not  essential  to  producing  a  high  quality  grid 
and  has  been  included  in  the  3-D  method  only  for  completeness.  The  second 
curvature  term,  S2,  is,  however,  critical  to  achieving  high  quality  grids 
around  the  leading  edge  of  the  delta  wing.  This  term  is  computed 
automatically  from  the  2-D  boundary  grids  and  then  interior  values  are 
interpolated  from  the  boundary  values. 

The  method  used  to  determine  the  values  of  SI  and  $2  in  the  interior  of 
the  3-L  grid  is  2-D  transfinite  interpolation,  more  correctly  known  as  blended 
bilinear  interpolation.  The  details  of  blended  bilinear  interpolation  are 
included  in  Appendix  C. 

The  3-D  elliptic  equations  are  solved  using  a  point  over-relaxation 
scheme.  Point  iteration  schemes  are  typically  slower  to  converge  than  are 
line  iteration  and  more  complex  schemes,  but  are  more  robust.  The  iteration 
process  requires  that  an  initial  guess  of  the  grid  point  locations  be  made 
prior  to  iterating.  The  favorable  characteristics  range  of  point  iteration 
schemes  permits  the  initial  guess  to  be  generated  by  simple  and  quick  methods. 
The  method  used  to  generate  the  initial  guess  is  blended  trilinear  interpola¬ 
tion.  The  details  of  blended  trilinear  interpolation  are  included  in  Appendix 
D. 

ALGEBRAIC  DOWNSTREAM  GRID 

The  grid  downstream  of  the  wing's  trailing  edge  is  generated  algebrai¬ 
cally.  A  2-D  grid  of  the  same  outer  shape  as  the  wing  trailing-edge  grid  is 
constructed  at  the  downstream  boundary,  consisting  of  concentric  ellipses, 
with  the  intermost  ellipse  being  a  slit.  line.  The  spacing  of  grid  points  on 
the  slit  line  and  the  spacing  between  the  concentric  ellipses  is  determined  by 


the  computer  code.  The  interior  of  the  downstream  grid  is  generated  by  simple 
linear  interpolation  between  the  wing  trai ling-edge  grid  and  the  downstream 
grid.  The  streamwise  spacing  of  the  interior  is  determined  by  the  hyperbolic 
tangent  stretching  routine.  The  grid  spacing  along  the  K-varying  lines  is 
continuous  across  the  wing's  trailing  edge. 

COMPUTER  PROGRAMS 

Two  computer  programs  perform  the  grid  generation  described  in  this 
paper.  The  first  program  generates  the  six  2-D  sides  (a  more  correct  term 
would  be  "surfaces")  of  the  3-D  grid  as  outlined  in  parts  1-5  of  this  section. 
Appendix  E  provides  directions  for  use  of  this  program.  The  second  program 
generates  the  complete  3-D  grid  and  requires  only  the  output  of  the  first 
program  as  input  data. 
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GRID  APPLICATIONS 

The  2-D  and  3-D  elliptic  grid  generation  codes  described  in  this  paper 
have  been  applied  to  several  geometries.  Two  configurations  of  interest  for 
which  Euler  solutions  have  been  obtained  are  a  cropped  delta  wing  and  a  blunt 
nose  missile  body.  The  cropped  delta  wing  configuration  is  depicted  in  Eigs. 
1-9.  Euler  solutions  for  this  configuration  at  Mach  numbers  of  0.4  and  0.85 
are  in  good  agreement  with  experimental  results.  Drag  polars  for  this 
configuration  (Figs.  10  and  11)  indicate  that  the  generated  grid  appears 
adequate  to  capture  the  major  flow  features,  including  what  appears  to  be 
vortex  breakdown.  The  predicted  drag  polar  has  been  shifted  to  reflect  zero 
lift  drag. 

The  delta  wing  grid  topology  was  applied  to  a  blunt  nose  elliptical  cross 
section  missile  body  (Ref.  5).  The  missile  nose  grid  and  a  slice  through  the 
3-D  space  is  shown  in  Fig.  12,  and  the  ?-D  missile  base  grid  (corresponding  to 
the  delta  wing's  trailing  edge)  is  shown  in  Fig.  13.  Again  the  3-D  grid 
generated  appears  adequate  to  resolve  the  major  flow  features  of  this  configu¬ 
ration  as  is  demonstrated  by  the  excellent  agreement  between  predicted  and 
experimental  normal  force  coefficients  at  Mach  2.0  for  angles  of  attack  up  to 
14.0  degrees  (Fig.  14). 
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SECTION  V 


CONCLUSIONS 

A  grid  topology  and  grid  generation  method  for  delta  and  cropped  delta 
wings  has  been  described  and  demonstrated.  The  grid  topology  has  been  adapted 
from  a  similar  topology  previously  used  for'  a  different  application.  The  grid 
generation  method  uses  a  blend  of  partial  differential  equations  and  algebraic 
interpolation.  The  resulting  grids  have  been  used  in  conjunction  with  an 
Euler  equation  solver  to  accurately  predict  the  flow  about  two  configurations. 
The  quality  of  the  prediction  is  extremely  comparable  to  results  produced  in  a 
similar  effort  of  advanced  grid  generation  and  flow  solution  for  a  delta  wing 
configuration  (Ref.  6). 
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Figure  5.  K  =  Constant  Grid  Layer  Showing  Rounded  Wingtip. 
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Figure  6c.  Plane-of-Symmetry  Grid  Showing  Root  Chord. 


Figurp  7.  Sketch  of  Trai  1  ing-Edge  Grid. 
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LIST  OF  SYMBOLS 

two-dimensional 

three-dimensional 

coefficient  of  drag,  referenced  to  wing  area 

coefficient  of  lift,  referenced  to  wing  area 

coefficient  of  normal  force,  referenced  to  base  area 

grid  indices 

maximum  1  grid  index 

maximum  J  grid  index 

maximum  K  grid  index 

grid  generation  source  term 

grid  generation  spacing  control  term 
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APPENDIX  R 


3-D  ELLIPTIC  EQUATION'S 
{Ref.  3) 
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BLENDED  BILINEAR  INTERPOLATION 
(Ref.  7  ) 

Purpose: 

Interpolate  function  F(s,t)  which  is  known  on  the  boundaries  of  a  2-D 
computational  space.  The  computational  space  has  dimension  which  vary  in  a 
direction  s  from  0  to  1  and  in  a  direction  t.  from  0  to  1. 

Assume  blended  bilinear  interpolants  of  the  form: 

Ps  =  (1  -  s)  F(0,t)  +  sF(l.t) 

Pt  =  (1  -  t)  F(s ,0)  +  tF(s,l). 

The  function  F(s,t)  is  interpolated  on  the  computational  space  by  the 


rule: 


F(s.t)  =  Ps  +  Pf  -  PsPt, 


F(s,t)  =  (1  -  s) 
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Grid  generation  is  actually  performed  on  a  computational  space  that 

varies  in  an  I  direction  from  1  to  I  and  in  a  J  direction  from  I  to  J 

max  max 

The  I  and  J  indices  are  related  to  s  and  t  by 
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BLENDED  TRILINEAP  INTERPOLATIONS 


(Ref.  7) 


Purpose: 


Interpolate  function  F(s,t,u)  which  is  known  on  the  boundaries  of  a  3-D 


computational  space.  The  computational  space  has  dimensions  which  vary  in  a 


direction  s  from  0  to  1,  in  a  direction  t  from  0  to  1  and  in  a  direction  u 


from  0  to  1. 


Assume  blended  tril inear  int.erpolants  of  the  form: 


Ps  =  (1  -  s)  F(0,t,u)  +  sF(l.t,u) 


Pt  =  (1  -  t)  F(s  ,0,u)  +  t.F(s,l,u) 


Pu  «  (1  -  u)  F(s,t,0)  '  uF(s.t.l). 


The  function  F(s,t,u)  is  interpolated  on  the  computational  space  by  *he 


rule: 


F(s,t,U)  -  P$  +  Pf  +  Pu  -  PsPf  -  PtPu  -  P/s  ♦  PsPtPt 


where. 


PPt  =  (1  -  s)  (1  -  t)  F(0,0,u) 

3  *  /  4  \/  .  \  r-  /  /-X  1  \ 


+  (1  -  s)  (  t)  F(0,l,u) 
+  (  s)  (1  -  t)  F( 1 ,0,u) 
+  (  s)  (  t)  F(1 ,1 ,u) 


PtP  =  0  -  t)  (1  -  u)  F(s ,0,0) 
f  U  +  (1  -  t)  (  u)  F(s  ,0,1) 

+  (  t)  (1  -  u)  F(s  ,1 ,0) 

+  (  t)  (  u)  F(s,l,l) 


P  P.  =  (1  -  u)  (1  -  s)  F(0,t ,0) 
U  5  Ml  -  u)  (  s)  F(  1  ,t ,0) 


+  {  u)  (1  -  s)  F(0,t,l) 

+  (  u)  (  s)  F( 1  ,t ,1) 


p  p*p 

S  t  l! 


S)  (1 
S)  (1 


s)  (1 
S)  (1 


t)  (1 
t)  ( 
t)  (1 
t)  ( 
t)  (1 
t)  ( 
t)  (1 
t)  ( 


u)  F(0,0,0) 
u)  F(0,0,1) 
u)  F(0, 1 ,0) 
u)  F(0,M) 
u)  F( 1 ,0,0) 
u)  F( 1,0,1) 
u)  F(l,l ,0) 
u)  F(l.l.l) 


The  relation  between  the  interpolation  variables  s,  t,  and  u  and  the  grid 


generation  variables  I,  J,  and  K  are: 


s  -  0  -  D  /  ('n»x  '  11 


‘  '  <J  -  1  <J„ax  ' 


u  =  (K  -  1)  /  (Kn;ax  -  1). 


APPENDIX  E 


2-D  Grid  Generation  Code  Input.  Directions 
Note:  all  data  is  in  free  format 

Note:  Cards  2,  4,  6,  10,  14,  19  and  22  are  data  set  header  cards  and  must 
be  input  as  indicated;  i.e.,  data  is  not  entered  on  those  cards. 


Card  # 


title 

MI 

MO 

MK 

#ITER 

SIGMA 

TOC 

mi 

mj 

mk 

£i  ter 

sigma 

toe 

mi  =  number  of  grid  pts  in  the  I  direction 

mj  =  number  of  grid  pts  in  the  J  direction 

mk  =  number  of  grid  pts  in  the  K  direction 

#iter  =  number  of  iteration 

sigma  =  spline  fit  coefficient 

0  -*■  cubic  spline  fit 
+  ®  -*■  straight  line  fit 

toe  =  percent  increase  of  input  airfoil  coordinates 


1  no  change  to  input  coordinates 


r  = 

relaxation 

success) 

factor. 

0  < 

r  < 

2  (1.1  has  been  used  with  some 

4. 

Tl 

T2  T3  T4 

T5  T6 

T7 

T8 

T9  T10  Til 

5. 

tl 

ro 

r-t 

CO 

<—► 

4^ 

t5  t6 

t7 

t8 

t9  tlO  til 

tl  =  percent  of  root  chord  arch  length  for  the  leading  grid  point 
spacing  on  the  root  chord. 

t2  =  percent  of  root  chord  arch  length  for  the  trailing  grid 
point  spacing  on  the  root  chord. 

t.3  -  spanwise  trai 1 ing-edge  grid  spacing  at  the  tip,  ratioed  to 
t.1;  t3  =  1.  will  cause  leading  spacing  to  be  constant  from 
root  to  tip. 

t4  =  spanwise  trai 1 ing-edge  grid  spacing  in  percent  semispan  at 
the  root  chord. 


tS  =  leading  edge  root  chord  stand-off  distance  of  J  =  2  layer, 
in  percent  root  chord. 

t.G  -  far-field  grid  spacing  along  J  line  extending  through  root 
chord  leading  edge,  in  percent  root  chord. 

t7  =  distance  from  leading  edge  root  chord  to  far  field  in 
percent  root  chord. 

tC  -  tip  spanwise  stand-off  distance  of  J  =  2  layer,  in  percent 
root  chord. 

t9  =  spanwise  grid  spacing  at  the  far  field  on  the  trailing  edge 
grid,  in  percent  root  chord. 

tlO  =  distance  spanwise  front  wingtip  trailing  edge  to  far  field, 
in  percent  root,  chord. 

til  =  vertical  distance  from  root  chord  trailing  edge  to 
far  field,  in  percent  root  chord. 

LEADING-EDGE:  SHAPE 

#lepts 

tflepts  =  number  of  points  defining  wing's  leading  edge  and 
tip. 

*1  *2  yn 

spanwise  coordinates  of  points  describing  the  wingtip  and 

leading  edge,  from  wingtip  trailing  edge  to  leading-edge  root 
chord. 

x,  x9  ...  xn 

id  n 

spanwise  coordinates  of  points  describing  the  wing  tip  and 

leading  edge,  from  wing  tip  trailing  edge  to  leading  edge  root 
chord. 

TRAILING-EDGE  SHAPE 
#tepts 

yl  y2  •••  yr. 

spanwise  coordinates  of  points  describing  the  wing  trailing  edge, 
from  root  to  tip. 

X1  x?  xn 

spanwise  coordinates  of  points  describing  thp  wing  trailing  edge, 
from  root  t.c  tip. 


UPPER  SURFACE  AIRFOIL  DATA 


15. 


ninsecu 


ninsecu  =  number  of  airfoil  sections  describing  the  wing  upper 
surface. 

16.  AIRFOIL  SECTION  IDENTIFIER 

17.  nptsec  xlesec  ylesec  xtesec  ytesec 

nptsec  =  number  of  data  points  describing  the  airfoil  section. 

xlesec,  ylesec  =  the  streamwise  and  spanwise  coordinate*"  of  the 

section  leading  edge. 

xtesec,  ytesec  =  the  streamwise  and  spanwise  coordinates  of  the 

section  trailing  edge. 

*1 

Z2 

Z3 


ZNPTSEC 

airfoil  section  coordinates  ncndimenional ized  by  local  chord, 
leading  edge  x  =  0. 
trailing  edge  x  =  1.  or  100. 

Repeat  cards  16,  17,  and  10  NINSECU  times  to  enter  all  upper  surface 
airfoil  sections. 

19.  WINGTIP  SHAPE 

20.  npttip  =  number  of  points  describing  the  true  wingtip  (zero 

thickness)  camber  shape. 


xnptt.ip  ^npttip 


10. 


NPTSEC 
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LOWER  SURFACE  AIRFUIl  L'A1 A 
ninsecl 

Repeat,  cards  1(> ,  17,  and  18  NINSECL  times  to  enter  al1  lower  surface 
airfoil  sections. 


PLANFORM  INPUTS 


(XTESEC'  YTESEC} 
AIRFOIL  SECTION  INPUTS 


Figure  El.  Input  Data  Definitions 
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